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NOMENCLATURE 


A 
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C 
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f 

F 
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h 
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m 

M 

P 

Pr 

%, 

R 

Re 
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constant,  Eqs.  (10)  and  (17) 
conversion  factor,  Eq.  (40) 
constant,  Eq.  (7) 

1  2 

friction  coefficient,  t  /■=■  p  u  ,  Eq.  (36) 

w  2  re  e  ^ 

conversion  factor,  Eq.  (40) 
conversion  factors,  Eq.  (40) 
velocity  function,  Eq.  (11) 
inviscid  flow  function,  Eq.  (8) 
conversion  factors,  Eqs.  (22)  and  (23) 
endialpy  function,  Eq.  (12) 
enthalpy,  Eqs.  (3)  and  (4) 
total  enthalpy 

exponent  of  power  law  shock,  Eq.  (7) 
normalized  lateral  flux,  Eq.  (39) 

Mach  number,  Eq.  (32) 

pressure 

Prandtl  number 

surface  heat  transfer,  Eq.  (35) 
inviscid  flow  function,  Eq.  (8) 

Reynolds  number,  Eq.  (33) 
function  defined  in  Eq.  (23) 
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S0.  Sv  S£,  Sj 
St 


t 

u 

v 
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x,  y,  z 
a 

Y 

6* 
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T) 

9 

fA 
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§ 

P 

o 

a 

T 

T 

W 

9 
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integral  functions,  Eqs.  (28),  (29),  (30),  and  (31) 
Stanton  number,  q/p  u  H  ,  Eq.  (37) 

W  C  C  6 

time 

velocity  component  in  x-direction 
velocity  component  in  y-direction 
transformed  lateral  velocity 
coordinate  axes 
(m-l)/m 

ratio  of  specific  heats 
displacement  thickness,  Eq.  (25) 
integral  thickness  defined  by  Eq.  (26) 
transformed  coordinate,  Eq.  (10) 
momentum  thickness,  Eq.  (27) 
viscosity 

kinematic  viscosity,  |±/p 

transformed  coordinate,  =  1  -  (x/xb),  Eq.  (10) 
density 

=  0,  two  dimensional;  =  1,  axisymmetric 

=  0,  plane  wave;  =  1,  axisymmetric;  =  2,  spherical 

transformed  coordinate,  Eq.  (10) 

surface  shear,  Eq.  (34) 

inviscid  flow  function,  Eq.  (8) 

stream  function,  Eqs.  (9)  and  (12) 
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outer  edge  of  boundary  layer 
initial  value 
wall  or  surface 
shock  front 
partial  derivative 

0(  )/8T),  alternative  notation  for  partial  derivation  of  T) 
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I.  INTRODUCTION 


The  solution  for  the  wall  boundary  layer  behind  a  shock  wave  moving 
with  nonuniform  velocity  is,  in  general,  complicated  because  of  the  addition 
of  the  temporal  dimension  to  the  problem.  However,  for  a  large  class  of 
flow  problems  of  practical  importance,  similarity  solutions  exist  so  that 
the  analysis  is  tractable.  Such  is  the  case  with  the  boundary-layer  flow 
behind  a  strong  shock  moving  with  power-law  velocity.  The  formulation  of 
]  the  problem  and  the  transformation  devised  to  reduce  the  independent  vari¬ 

ables  from  three  (x,  y,t)  of  the  original  problem  to  two  normalized  spatial 
variables  (S,T|)  in  the  transformed  frame  have  been  detailed  in  Ref.  1. 
Numerical  solutions  were  also  obtained  in  Ref.  1  by  series  expansions  of 
both  the  inviscid  and  boundary-layer  flow  in  powers  of  the  normalized  dis¬ 
tance  from  the  shock  wave,  §.  Two  terms  in  each  expansion  were  used,  and 
the  results  are  limited  to  the  region  directly  behind  the  shock  wave  (§2«  1). 
The  effort  involved  in  extending  the  series  to  higher  order  terms  increases 
very  rapidly.  Chen  and  Chang,  ^  apparently  unaware  of  the  existence  of 
Ref.  1,  reformulated  the  problem  using  essentially  the  same  approach,  but 
they  carried  the  expansion  to  the  next  term.  Unfortunately,  some  of  their 
coefficients  appear  to  be  in  error. 

Numerical  finite  difference  methods  are  available  for  evaluating 

3 

boundary-layer  flows  in  two  coordinate  variables.  We  have,  therefore, 
reexamined  the  shock-induced  boundary-layer  problem  in  the  light  of  a 
direct  application  of  finite  difference  methods  to  the  transformed  boundary- 
layer  equations.  In  this  report  similarity  (§,11)  solutions  are  obtained  for 
the  laminar  boundary  layer  behind  a  power-law  shock  associated  with  a 
'  blast  wave.  A  finite  difference  method  based  upon  Blottner’s  numerical 

t  scheme  is  used.  The  results  are  valid,  at  all  times,  in  the  entire  flow 

| 

I  region  between  the  shock  front  and  the  immediate  vicinity  of  the  blast-wave 


origin  provided  the  boundary  layer  remains  laminar.  The  method  of 
analysis  is  described  in  Section  II.  Results  and  discussions  are  given  in 
Section  III.  The  relation  between  blast-wave  strength  and  the  energy  of  the 
disturbed  flow  is  discussed  in  Appendix  A.  The  extent  of  the  laminar 
boundary  layer,  behind  blast  waves,  is  discussed  in  Appendix  B. 


II.  ANALYSIS 


Similarity  properties  of  the  inviscid  flow  associated  with  strong 
power-law  shocks  and  the  corresponding  similarity  transformation  devised 
for  the  shock-induced  laminar  boundary  layer  have  already  been  detailed  in 
Ref.  1.  For  completeness,  the  governing  equations  and  basic  assumptions 
are  briefly  repeated  here. 

With  the  assumption  that  the  fluid  Prandtl  n timber  and  specific  heats 
are  constant  and  that  the  viscosity  is  proportional  to  temperature,  the 
unsteady,  compressible  laminar  boundary-layer  equations  for  a  perfect  gas 
are: 

Continuity 


9£  +  _L 

at  a 

x 


Momentum 


Du 

P  Dt  ~ 


Energy 


Dh  Dpe 
p  Dt  '  Dt 


+  3(£yl  =  o 

dx  T  ay 


State 


p  =  [(Y  *  1)/Y]ph 


(1) 


(2) 


(3) 


(4) 
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where 


Dt  '8t  +  8x  +  V  ay 

Symbols  are  defined  in  the  Nomenclature  section,  and  the  orientation  of 
coordinates  is  indicated  in  Fig.  1.  Both  two-dimensional  and  axisymmetric 
boundary  layers  corresponding  to  o  -  0  and  1,  respectively,  are  included. 
The  choice  of  a  is  determined  by  the  geometry  of  the  propagating  shock  and 
the  surface  on  which  the  boundary  layer  develops  (Fig.  2).  The  boundary 
conditions  are: 


at  y  =  0  u  (x,0,t)  =  0 

v  (x,  0,t)  =  v^  (x,  t)  (5) 

hw  (x,  0,t)  =  h^  (x,t) 

at  y  =  ®  u  (x,®,t)  =  ug  (x,t) 

h  (x,»,t)  =  hg  (x,t)  (6) 

We  confined  our  attention  to  a  strong  shock  that  moves  according  to 
the  power  law  (Appendix  A) 


x  =  Ct 
s 


m 


u  = 

s 


dx 
_ 8 

dt 


Cmt 


m- 1 


(7) 


The  disturbed  flow  for  constant  specific  heats  is  similar  at  successive 

instants  in  time.'*  If  the  independent  variable  x  is  replaced  by  5  =  1  -  (x/x  ), 

s 

the  disturbed  flow  can  be  expressed  as 


(8) 


M  1  BWW' 


wm 


p./p. 


p./p.  -  R<s> 


"e  =  “,*(5). 


_x_  5a  =_*_  .2  HSi 

Y  -  1  Pe  Y  -1  8  R(§) 


The  functions  F,  R,  and  cp  depend  on  the  geometry  of  the  shock  wave  as  well 
as  the  parameters  y>  m>  and  the  coordinate  §.  Analytical  and  numerical 
solutions  for  F,  R,  and  (p  are  given  in  Ref.  4  for  plane,  cylindrical,  and 
spherical  shocks  (corresponding  to  o'  =  0,  1,  2  respectively)  and  y  =  1,4. 
When  m  =  2/(3  +  o),  the  flows  correspond  to  constant  energy  blast  waves. 
These  cases  are  discussed  in  Appendix  A. 

Four  combinations  of  shock  and  surface  orientation  are  of  practical 
interest.  These  are  shown  schematically  in  Fig.  2. 

Equation  (1)  is  satisfied  by  a  scalar  stream  function  ^  such  that 


p« 

fi , 

< 

II 

1 

k 

k  +  jl  I, 

,p/-V  A] 

o 

ay 

a 

dx  at  I 

px 

px 

L  V 

*/o  *  /J 

The  following  independent  variables  are  introduced; 


g  =  l-  x/Ctm 

,  (t) d7 

[At2m(CT+1)  -  la  -x/ctm)]iyz 

T  =  t 

where  A  is  a  constant  to  be  defined  later.  Next,  the  following  form  is 
assumed  for  the  dependent  variables 
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(11) 


*  =  u  [AT2m(0+  1}  ■  15]1/2f (S,  T)) 

© 

h/h  =  p  /p  =  g(g,  1))  (12) 

©  © 

which,  by  Eq.  (9),  gives  f^  =  u/ue«  For  the  present,  the  surface  is  assumed 
impermeable.  The  boundary  conditions  on  f(?,  11)  are  then 

f(§»  0)  =  fjj  (5,  0)  =  0,  £^(5,  •)  =  1  (13) 

Wall  temperature  is  assumed  negligible  compared  with  the  free- stream 
temperature.  The  boundary  conditions  on  g  (?,  11)  are  then 

g  (?.  0)  =  0,  g(5,  »)  =  1  (14) 

When  Eqs.  (10),  (11),  and  (12)  are  introduced  into  Eqs.  (1),  (2),  and 
(3),  they  become 

(1  -  ;>2°<F/Fo)fwt 

■  +*fj  -4  (2°  + 

(1  -  5)2°^  Kjjjj  t^T|| 

+  (11  -  <pf)  g-jj  =  2?  j[fcp5  +  cpf  ?  (2o  +  cr)  lljg^ 


(16) 


where  a  =  (m-l)/m.  In  deriving  Eqs.  (15)  and  (16)  the  constant  A  was 
chosen  as  follows  to  simplify  the  coefficients  of  f^^  and  g^ 

A  =  2mF  C2(a+1)v  p  /p  (17) 

O  00  00  00 

The  explicit  dependence  of  the  flow  on  the  independent  variable  t  has  now 
been  eliminated. 

The  development  up  to  here  is  the  same  as  that  presented  in  Ref.  1. 
Equations  (15)  and  (16)  are  two  simultaneous  equations  for  f  and  g  in  terms 
of  the  two  independent  transformed  variables  §  and  11.  Functionally,  they 
are  analogous  to  the  steady- state  boundary  layer  equations  successfully 
treated  by  Blottner.  Blottner's  derivation,  however,  made  use  of  a  trans¬ 
formed  lateral  velocity  component  that  reduced  the  order  of  the  differential 
equation  for  f  by  one.  The  lateral  velocity  component  is,  in  turn,  defined  by 
an  additional  differential  equation  arising  from  the  continuity  equation. 

The  transformed  lateral  velocity  component  is  defined  by 


V  =  -2  -a 
<P 


which  is  motivated  by  the  expanded  form  for  ipx-  The  relation  between  V 
and  v  is  found  with  the  use  of  Eq.  (9).  Eqs.  (15),  (16),  and  (18)  now  become 

a?  -»]!£-  [*v  *  ti  |25  (<.+!)  +  ijj  f  t\ 

+  25^  +  11-?)^  -*{f'  -if  ?']  •  f  f(pr 0 


25  <1-5-5*')  |f  -  [*v  ♦  H  j25(o  +f)  +  lj]g„ 

-  ^  «P2  #■  +  25  j(l-5-»f')  f1] 


..  ..2c  FI  _  A 

(1-5)  F  pr  -  0 

o 


(20) 


T^25^*1)  f'  +  2E!r=0  (2,) 

Note  that  i'  (s3f/3Tl)  is  now  the  lowest  order  of  the  function  f  appearing  in 
the  above  equations,  which,  termwise,  are  now  directly  comparable  to 
Blottner's  equations  of  Ref.  3. 

A  computer  program  has  been  written  based  on  the  numerical  scheme 
outlined  in  Ref.  3,  and  Eqs.  (19),  (20),  and  (21 )  have  been  integrated  for  the 
four  cases  of  practical  interest  shown  in  Fig.  2.  The  results  are  presented 
in  Section  III.  Quantities  of  interest  are  noted  below. 

Physical  coordinates  are  obtained  from  [e.g.  ,  Eqs.  (10)  and  (17)] 

x  =  xg  (1  -  §) 

y=F1-F2.S  (22) 


where 
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F,<"  =  [Fo  jT 


F2<!>  --»*»£ - 
1  (l-fi)R 


5(5,11)  =/‘11g  dTl 


Boundary-layer  integral  lengths  are  expressible  as: 


ye  =FrF2*S0  (5) 


-f  (!  *p^) dy "  FrF2'si(5) 


dy  =FrF2.52(5) 


Fi-F2’S3(5) 


where 


r^e  „ 
S0  =J  g  dH 
J0 

r\ 

\  f  <■-*'! 

*'0 


<•#**»*-  ■-< 


-  .  J~»  •  *,  *«t*»*i.  • 


(g-1)  dl) 


(30) 


*'0 

=  /  *  f'  <i-f') 


Here  11  is  the  value  of  1)  at  the  outer  edge  of  the  boundary  layer,  where  the 
6 

edge  boundary  conditions  (y  =  °°)  are  satisfied  in  the  present  numerical  code. 
The  choice  for  11  is  discussed  later.  The  quantity  yg  is  the  value  of  y 
corresponding  to  1)  and  is  a  rough  measure  of  boundary  layer  thickness. 

jj,  C  jjc  >J< 

The  quantity  6  is  the  familiar  displacement  thickness.  The  use  of  6  will 
be  discussed  later;  0  is  the  momentum  thickness.  Note  that  these  integral 
lengths  are  dominated  by  the  factor  F^,  which  is  singular  at  (  =  i. 

Local  Mach  number  M  and  Reynolds  number  Re  can  be  expressed  in 
the  form 


(ypjpe) 


p  u  (x  -  X 
Ke  e  s 


(yF/R) 


The  shear  and  heat  transfer  at  the  wall  can  be  found  from 


t  =  p. 
w  w 


1/2 

Mp.Vn<fc0'  <“> 


(35) 


“w  "  -  kw  (” )w  =  ■  Pr  [2PoP;».  x  J  [^1^]  P«h«*t|  < '•  °> 


Normalized  expressions  for  the  above  are 


°f  <*«>1/2  (Re)1/2 -(w-)in  (1-5)<I  fW5,0)  (361 


2  re  e 


St  (Re) 


1/2 _ qw _  lnA/2 


w 

p  u  H 
^e  e  e 


(Re) 


9(1  +  rTrrrrl 


2  yF/R 


) 


(37) 


where  Cf  and  St  represent  the  friction  coefficient  and  Stanton  number, 
respectively.  These  are  functions  only  of  ?. 

Integration  of  the  continuity  equation  [Eq.  (1)]  across  the  boundary 
layer  provides  an  expression  for  net  mass  flux  in  the  lateral  direction, 
namely 


«<">.  -  ■*.«.** 


(38) 


which  can  be  expressed  in  the  form 


^s(25)1/2  [(pv)e  -  (pv)w]/B 

r  dS2  dSll 

=  c2  [D1S2  +  V  -df  -  D3S1  -  d4  -3fJ 


(39) 
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where 


B  =  p  [V  p.  u^/p  x  I 
L  °  ®  »  F®  sj 


1/2 


C2 

=  1/(1- 5)° 

Di 

=  ‘t25(f 

D2 

=  25  (1-5) 

D3 

=  25<P?+  <P 

°4 

=  25<P 

(40) 


The  normalized  lateral  mass  flux  is  a  function  only  of  §.  The 
present  computing  program  provides  Sj  and  S2  versus  5  in  appropriate 
tables  such  that  the  derivatives  dSj/d5>  dS2/d§  and  the  corresponding  values 
oi~4(  can  be  subsequently  calculated. 
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SHOCK 


X  xs  / *s 

i  z 


(a)  PLANE  MOVING  SHOCK  (b)  CYLINDRICAL  MOVING  SHOCK 
(cr  =  0,  <T  =  0)  (cF=1,  (T  =  0) 


Y 


Ic)  CYLINDRICAL  MOVING  SHOCK  (d)  SPHERICAL  MOVING  SHOCK 
Icr  =  1,  <7  =  1)  (<r  =  2,  cr=  1) 


Fig.  2.  Moving  shock-plane  wall  orientation 
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HI.  RESULTS  AND  DISCUSSIONS 


Inviscid  and  viscid  flow-field  solutions  have  been  obtained;  the  results 

are  discussed  herein.  The  inviscid  flow  functions  F,  R,  and  cp  associated 

with  blast  waves  are  discussed  in  Ref.  4.  Both  analytical  expressions  and 

numerical  tables  for  the  cases  c  =  0,  1,  2,  and  y  =  1.4  are  given.  For  the 

present  applications,  however,  it  has  been  found  necessary  to  refine  the 

tabular  results  to  minimize  the  interpolation  errors  in  coefficients  appearing 

in  the  boundary-layer  equations.  To  this  end,  it  was  found  more  convenient 

to  recalculate  the  inviscid  flow  functions  and  their  derivatives  to  the  degree 

of  accuracy  desired  by  starting  from  the  basic  equations  governing  a  power- 

law  shock,  such  as  those  given  in  Ref.  5.  The  results  for  y  =  i.4  are  given 

—  ❖ 

here  as  Tables  1,  2,  and  3,  for  the  cases  o  =  0,  1,  and  2,  respectively. 

They  are  also  presented  as  Figs.  3,  4,  and  5.  The  variation  of  Mach  num¬ 
ber  M  with  5  for  the  three  inviscid  blast-wave  flows  are  plotted  as  Fig.  6. 
The  variation  of  Re/(p  x  /p.  u  )  with  5  for  the  same  cases  are  plotted  as 
Fig.  7. 

Boundary-layer  solutions  have  been  obtained  for  y  =  i.4,  Pr  =  0.  72 
for  the  four  cases  noted  in  Fig.  2.  In  each,  the  integration  has  been  suc¬ 
cessfully  carried  out  from  the  shock  to  the  immediate  vicinity  of  the  blast 
origin.  At  the  blast  origin  (§  =  1.0),  both  the  inviscid  flow  and  the  boundary- 
layer  equations  are  singular.  Results  have  been  obtained  for  all  cases  to 
well  beyond  ?  =  0.  9.  Divergence  did  set  in,  however,  as  the  integration 
approached  5=1.  The  refinement  of  the  inviscid-flow  functions  discussed 
earlier  contributed  noticeably  to  extending  the  range  of  integration  closer  to 
the  singularity.  The  numerical  techniques  and  the  results  are  discussed  in 
the  following  paragraphs. 


$ 

The  authors  are  indebted  to  Heather  Bagwell  of  The  Aerospace  Corporation 
for  providing  these  results. 
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The  boundary-layer  equations  were  programmed  using  Blottner's 
3 

numerical  technique.  At  each  streamwise  station  5  the  present  computing 

program  allowed  up  to  100  mesh  points  distributed,  either  uniformly  or 

nonuniformly,  in  the  interval  0  <  T)  £  7]  ,  where,  as  previously  noted,  7) 

denotes  the  value  of  1]  at  which  the  free  stream  boundary  conditions  (y  =  <*>) 

are  applied.  For  the  results  herewith  presented,  98  mesh  points  were  used. 

A  graduated  grid  size  was  used  so  that  the  mesh  network  was  finest  near 

the  surface  where  the  definition  of  the  profiles  is  of  the  greatest  interest. 

The  calculation  is  most  accurate  when  the  value  of  7)  at  each  streamwise 

e 

station  is  just  sufficient  for  the  profiles  to  approach  the  outer  edge 
asymptotically  within  an  acceptable  tolerance. 

The  boundary-layer  solution  at  5  =  0  was  obtained  by  iteration,  using 
reasonable  profiles  for  f^  and  g  as  first  estimates.  For  5=0,  Eqs.  (19) 
and  (20)  reduce  to  the  zero-order  equations  of  Ref.  1;  the  two  results  are, 
therefore,  directly  comparable.  As  has  been  noted  in  Ref.  1,  they  do  not 
depend  on  a  or  a.  The  wall  gradients  of  the  present  calculation  are  listed 
in  Table  4  and  are  shown  to  compare  very  favorably  with  those  of  Ref.  1. 

A  first  estimate  of  the  boundary-layer  solution  for  ?  >  0  was  obtained 
by  satisfying  outer-edge  boundary  conditions  at  1)  =6  for  all  values  of  5. 

The  choice  T)  =  6  was  based  on  the  5=0  solution.  The  resulting  velocity 
and  enthalpy  profiles  at  successive  values  of  5  are  given  in  Fig.  8  in  the 
form  of  composite  profile  summary  plots.  Note  that  the  velocity  profiles 
develop  overshoot  with  increase  in  |.  It  is  also  seen  from  Fig.  8  that,  in 
terms  of  T|,  the  boundary-layer  thickness  decreases  significantly  with 
increase  in  5.  That  is,  the  physical  edge  of  the  boundary  layer  occurs  at 
significantly  smaller  values  of  T1  as  5  increases.  (This  is  a  consequence 
of  the  definition  of  T|.  In  the  physical  variable,  y,  the  boundary- layer  thick¬ 
ness  increases  with  increase  in  5-  )  In  order  to  maintain  computational 
accuracy,  it  is  desirable  to  decrease  1)  with  increase  in  5*  The  present 
approach  was  to  let  the  second  choice  for  7)  correspond  to  the  values  of  7), 


at  which  u/ue  s  =  1  ±  0.  0005  in  the  preceding  T|^  =  6  solution.  The  +  and 
-  signs  refer  to  velocity  profiles  with  and  without  overshoot.  The  resulting 
values  of  H  are  given  in  Fig.  9.  The  boundary-layer  profiles  were 
recalculated  using  the  new  -distribution.  These  calculations  not  only 
provided  a  more  detailed  description  of  the  profiles  by  utilizing  more  fully 
the  effective  Tl-thickness  of  the  boundary,  but  also  improved  the  stability  of 
the  numerical  process  at  large  §  where  the  TJ-thickness  can  be  as  much  as 
two  orders  of  magnitude  smaller  than  that  at  ?  =0.  It  was  found,  however, 

that  the  results  for  wall  shear  and  heat  transfer  and  for  the  boundary  layer 

*  ## 

thicknesses  6,6  ,  and  9  from  the  second  calculation  did  not  differ  signi¬ 

ficantly  from  the  first  T]  =  6  calculation  (i.  e.  ,  were  within  computational 
accuracy).  Hence  further  iterations  of  T]  were  not  needed. 

The  results  from  the  variable  1)  solution  are  given  in  Figs.  10  to  13 
and  in  Tables  4  and  5.  Sample  velocity  and  enthalpy  profiles  are  given  in 
Fig.  10.  These  figures  provide  somewhat  more  detail  than  Fig.  8.  Note 
that  the  velocity  overshoot  is  least  severe  for  the  case  o  =  a  =  0  and  is  most 
severe  for  the  case  a  =  1,  a  =  0.  A  maximum  overshoot  velocity  of  about 
u/u^  =  1. 09  is  achieved  at  ?  =  0.  5  in  the  latter  case.  Velocity  and  enthalpy 
gradients,  evaluated  at  the  wall,  are  given  in  Fig.  11  and  Table  5.  The 
derivative  of  f^^  (§,  0)  and  g^  (5,0)  with  5  at  %' =  0  has  been  obtained  numeri¬ 
cally  in  this  study,  and  a  comparison  is  made  in  Table  4  with  the  correspond¬ 
ing  values  from  Ref.  1.  Agreement  is  to  within  about  three  significant 
figures,  which  is  about  the  accuracy  of  some  of  the  present  computer  code 
subroutines.  Wall  shear  and  wall  heat  transfer  coefficients,  referenced  to 
local  flow  properties  [Eqs.  (36)  and  (37)],  are  given  in  Fig.  10.  It  is  seen 
that  the  shear  and  heat  transfer  coefficients  increase  by  about  an  order  of 
magnitude,  as  5  increases  from  0  to  about  0.  9.  The  results  for  normalized 
lateral  mass  flux  and  boundary-layer  lengths  are  given  in  Table  5  and 
Fig.  13.  These  quantities  have  been  normalized  to  remain  of  order  one  as 

5  varies.  The  value  of  S_  in  Table  5  is  based  on  the  variable  T)  solution 

0  e 

(Fig.  9).  Hence  the  corresponding  value  yg  provides  a  rough  estimate  of 


local  boundary- layer  thickness.  Equations  (22)  to  (31)  indicate  that  in 
physical  variables  the  characteristic  boundary-layer  lengths  (e.  g.  , 

★  Q 

ye#6  ,  .  .  .)  are  proportional  to  1/(1-?)  R.  Because  R  —  0  as  ?  -—  1 ,  these 
lengths  tend  to  become  infinite  as  ? -*•  1,  for  a  =  0  as  well  as  o  =  1  cases. 

Because  both  the  inviscid  and  the  boundary-layer  flow  fields  near  the 
blast  origin  are  more  or  less  idealized,  the  solution  in  this  regime  is  more 
of  analytical,  than  practical,  interest  in  that  it  serves  as  a  test  for  the 
functional  self-consistency  of  the  equations  and  the  method  of  calculation 
used. 
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Table  1.  Inviscid  flow  functions, a  =  0,  y  =  1.4  (Continued) 
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Table  2.  Inviscid  flow  functions,  d 
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Table  2.  Inviscid  flow  functioned  =  1,  y  =  1.4  (Continued) 
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Table  3.  Inviscid  flow  functions,  0  =2,  Y  =  (Continued) 
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Table  4.  Comparison  with  wall  gradients  of  Reference  1 


V  =  1.  4 

Pr  =  0.  72 

f"  =  (0,0) 

g'  (0.0) 

Present 

Ref.  1 

Present 

Ref.  1 

0. 66198 

0.66141 

0.89864 

0.  89693 

■ 

1 

a 

[9f  #/35]g_0 

[agV95]?=0 
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Present 

Ref.  1 

Present 

Ref.  1 

0 

1 

1 

2 

5 

-1/2 

-1 

-1 

-3/2 

1.  785 

2.  597 

3.  775 

4.  576 

1.  787 

2.  601 

3.  788 

4.  601 

-1.435 

-3.025 

-1. 814 

-3.405 

-1. 437 

-3.040 

-1.816 

-3.419 

7*  af; 

1  fw  (0,0) 


■  g'  (0,0) 
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Fig.  8.  Composite  profile  summary,  Case  B 


Fig.  8.  Composite  profile  summary.  Case 
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Fig.  10A-1.  Detailed  boundary-layer  profiles,  C  =0.10,  Case  A 
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Fig.  10A-3.  Detailed  boundary-la^ 
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10A-4.  Detailed  boundary -layer  profiles,  C  =  0.40,  Case  A 


10A-6.  Detailed  boundary-layer  profiles,  C  =  0.  60,  Case  A 


iOA- 


Fig.  10A-9.  Detailed  boundary-layer  profiles,  C  =0.90,  Case  A 


10B-2.  Detailed  boundary-layer  profiles,  Q  =  0.20,  Case  B 


er  profiles,  C  =  0.30,  Case  B 
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CASE  C 
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10C-3.  Detailed  boundary-layer  profiles,  £  =  0.30,  Case  C 


CASE  C 


Fig.  iOC-5.  Detailed  boundary-layer  profiles,  C  =0.50,  Case 


dary-layer  profiles,  Q  =  0.  60,  Case  C 


iOC-7.  Detailed  boundary-layer  profiles,  Q  =  0.  70,  Case 


Fig.  10D- 


Fig.  10D- 


Fig.  10D-4.  Detailed  boundary-layer  profiles,  £  =0.40,  Case 
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10D-5.  Detailed  boundary -layer  profiles,  Q  =  0.  50,  Case 
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Fig.  11.  Wall  gradient  of  f'and  g,  Case  B 
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Fig.  11.  Wall  gradient  of  f'  and  g,  Case  D 
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IV.  CONCLUDING  REMARKS 


The  results  of  this  study  provide  the  first  exact  solutions  of  the 

laminar  boundary  layer  behind  power -law  shocks  that  cover  the  entire 

terrain  swept  by  a  blast  wave,  with  the  exception  of  a  small  area  near  the 

center  of  the  blast  where  the  inviscid  and  viscid  equations  are  singular. 

Upon  increase  in  g,  the  results  depart  significantly  from  the  previous 

solutions  for  g^«  1.  In  particular,  in  all  four  cases.  Case  A  through 

Case  D,  the  normalized  heat  transfer  showed  a  minimum  before  g  =0.5. 

For  the  axisymmetric  boundary  layers  (®  =  1,  Case  C  and  D)  both  wall 

shear  and  heat  transfer  increase  at  much  faster  rates  with  increase  in  g 

2 

than  were  indicated  by  the  solutions  for  g  «  1. 

Boundary-layer  transition  is  discussed  in  Appendix  B.  The  present 

results  are  applicable  only  in  that  portion  of  the  flow  field  where  the 

boundary  layer  remains  laminar.  For  strong  shocks  in  air  or  argon,  the 

-3  -2 

boundary  layer  is  laminar  for  p  x  S  0  (10  -10  )  atm-ft.  For  larger 

values  of  p  x  ,  the  fraction  of  the  disturbed  field  that  remains  laminar 
3  -2 

is  g.  =  0(10"  - 10  )/p  x  ,  where  p  x  is  in  units  of  atm-ft.  Therefore, 

the  results  presented  here  are  of  primary  interest  for  low  pressure  test 
facilities  and,  possibly,  for  studies  of  the  interaction  of  pulsed  laser 
induced  blast  waves  with  ambient  surfaces. 
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APPENDIX  A 


BLAST- WAVE  STRENGTH 


The  variation  of  blast-wave  radius  x  with  time  t  is  discussed  herein. 

8 

For  blast  waves  with  constant  energy  E  the  variation  of  shock  radius  with 
time  can  be  expressed^ 


where  a  is  a  nondimensional  constant  that  depends  on  y  and  o.  Numerical 

_  4  5 

estimates  for  a  are  obtained  from  ' 

"aA(sf 5°d5  (A-2 

where  A  =  2,  2ir,  and  4ir  for  o  =  0,  1,  and  2,  respectively.  The  choice 
A  =  2  for  o'  =  0  indicates  that  the  latter  corresponds  to  a  symmetric  blast 
(i.  e.  ,  -1  <  5  <  1).  The  units  of  E  in  Eq.  (A-i)  are  energy/area,  energy/ 
length,  and  energy  for  o'  =  0,  i,  and  2,  respectively.  Numerical  estimates 
for  a  are  provided  in  Table  A-i. 
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from  Ref.  5_  (o  =  0,  ij  and  H.  Bagwell  (o  =  2) 
Values  for  cr  =  0  correspond  to  symmetric 
blasts,  -Is  5  s  1. 


APPENDIX  B 


BOUNDARY-LAYER  TRANSITION 


The  boundary  layer  is  laminar  directly  behind  the  shock,  but  it  may 
undergo  transition  to  a  turbulent  boundary  layer  at  some  distance  behind  the 
shock.  The  location  of  the  transition  point  defines  the  extent  of  the  laminar 
boundary  layer  and,  therefore,  the  region  of  validity  of  the  present  theory. 
The  latter  is  discussed  herein. 

Assume  that  local  Reynolds  number  Re  can  be  characterized  by  flow 
conditions  directly  behind  the  shock.  An  appropriate  Reynolds  number 

(based  on  distance  a  freestream  particle  has  traveled  relative  to  the  wall) 

.  6 

is 


l 


(B  - 1 ) 


where  R(0)  is  the  density  ratio  across  the  shock  and  v  =  |a  /  p  is  evaluated 

-2  e  e  eu) 

at  5  =  0.  For  an  ideal  gas,  a  strong  shock  (Ms  «  1)  and  n  —  T  ,  Eq.  (B-i) 
becomes 


(B-2) 
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For  most  gases,  u>  =  1/2.  Thus  the  dependence  of  Eq.  (B-2)  on  shock  Mach 

number  M  is  weak.  Taking  to  =  1/2  in  Eq.  (B-2)  yields 
8 


Re _ 4(y+l)  Pg>a<» 

^  Xs  [2y(V“  1 ) 5] 


(B-3) 


Let  and  Re^  denote  the  transition  location  and  the  transition  Reynolds 
number,  respectively.  For  air  and  argon,  Eq.  (B-3)  can  be  expressed  in 
the  form 

=  2.  5  X  10  ^  atm-ft  (air)  (B-4a) 

=  8.  4  X  10  ^  atm-ft  (argon)  (B-4b) 


p  x  | 


lot 

* lRe‘  T» 


where  x  =  x  -  xt  denotes  the  distance  behind  the  Bhock  at  which  transi- 

tion  occurs  and  is  a  direct  measure  of  the  fraction  of  the  flow  field  which 

is  laminar.  The  bracketed  term  on  the  left-hand  side  of  Eq.  (B-4)  is  of 

order  one.  Hence,  the  boundary  layer  is  laminar  in  the  entire  disturbed  flow 

-3  -2 

region  [i.  e.  ,  =  0(1)]  for  p^Xg  s  0  (10  -  10  )  atm-ft.  The  latter  regime 

occurs  primarily  in  test  facilities  (e.  g.  ,  blast-wave-driven  low  pressure 

shock  tubes)  but  may  also  be  of  interest  in  connection  with  the  study  of  pulsed 

laser  induced  blast  waves.  With  increase  in  p  x  above  values  of  the  order 
-2  “  s 

of  10  atm-ft,  the  fraction  of  the  disturbed  flow  that  is  laminar  is  reduced 

(i.  e.  ,  S«.~—0).  For  ?  <  1  and  a  given  ambient  gas,  x  -  x  depends  only  on 
t  t  8  t 

p^  and  is  independent  of  blast  energy  and  shock  radius. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos¬ 
pheric  optics,  cliemi cal  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro -optics;  communication  sciences,  applied  electronics,  semi¬ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia¬ 
tion  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  space  systems. 
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